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The effect of small, periodic phase errors upon transmission between 
two coaxial, circularly symmetric apertures is considered when the aperture 
phase distributions are confocal and the amplitude distributions are 
gaussian. The results are applicable to loss calculations in beam wave- 
guide systems with imperfect lenses. When the periods of the phase errors 
are less than one-half the aperture radii, the total loss is approximately 
|(/3i + /3i), where @i, /3 2 are the peak phase errors (in radians) on the 
apertures. Phase errors with periods greater than the aperture diameters 
are found to cause comparatively little transmission loss. 

I. INTRODUCTION 

The use of beam waveguide 1 systems for the transmission of informa- 
tion, 2 or for the transmission of power, 3 necessitates the design of lenses 
(or cylindrical reflectors 4 ) as focusing elements. In the design of these 
elements, it is desirable to estimate the degradation in performance 
caused by surface profile errors. Such degradation results in trans- 
mission loss and, in a communications system, will contribute to 
interference. Typically, the profile errors are associated with machining 
operations and, for lenses with circular symmetry, these errors are 
frequently circularly symmetric. The principal effect of the errors 
is to impart small, circularly symmetric phase perturbations to the 
field distribution adjacent to the lenses. The purpose of this paper 
is to calculate the reduction in transmission, caused by phase errors of 
this type, in a simple system comprising two coaxial, circular apertures 
as shown in Fig. 1. The field distributions on the apertures may repre- 
sent the fields in the aperture planes of two antennas or the fields on 
adjacent lenses in a beam waveguide system. 

In the absence of phase errors the transmission between coaxial 
apertures has been extensively studied by Kay, 5 Borgiotti, 6 Heurtley, 7 
and others, with the principal objective of determining that aperture 
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Fig. 1 — Coaxial, circular apertures with confocal field distributions S,'(r,') and 
ES(rJ). 



field distribution which maximizes the transmission. When the aperture 
separation is greater than some minimum, it has been found that this 
field distribution corresponds to that of the lowest order mode in an 
open, confocal resonator. 8 The phase distribution appropriate to this 
mode is obtained when the phase fronts on the apertures are confocal, 
i.e., are spherical with the center of curvature at the center of the 
other aperture. The appropriate amplitude distribution is well approxi- 
mated 9 by a gaussian curve. For this (optimum) distribution, the 
transmission between the apertures can attain surprisingly high values 
even when the apertures are separated by many aperture diameters 
(see Refs. 5, 6, and 7). Although the effect of periodic and random 
phase errors upon antenna gain and side lobe level has been investi- 
gated by several authors, 1011 little information appears to be available 
regarding transmission between two apertures when each has phase 
errors. In the case of transmission between two reflector antennas, 
Chu 12 has obtained an upper bound for the loss resulting from those 
phase errors produced by displacement of the feeds from the reflector 
foci. Yoneyama and Nishida 13 have considered transmission through 
a two-dimensional, confocal beam waveguide consisting of lenses with 
random phase errors. We compare their results to those of the present 
study later in this paper. 

In the following section the total transmission loss in a confocal 
system, with small phase errors on apertures with arbitrary ampli- 
tude distributions, is expressed in terms of the losses associated with 
each aperture when the other is free from phase errors. In the next 
section explicit expressions are derived, for two cases of practical 
interest, when the phase errors are sinusoidal and when the aperture 
amplitude distributions are gaussian. These expressions are then 
discussed and compared with results from the literature. 

784 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1975 



II. TRANSMISSION BETWEEN CONFOCAL APERTURES WITH 
PHASE ERRORS 

Consider the circular apertures A h A 2 of radius oi, a 2 which are 
separated by a distance d > ai, a 2 , as in Fig. 1. It is assumed that the 
apertures are focused at each other such that the tangential electric 
fields in the apertures, when each is transmitting in the absence of the 
other, have the quadratic phase variation 

SKu) = EM) exp (j |£ ) , i = 1, 2, (1) 

where the r\ are radial coordinates in the apertures. In the absence of 
phase errors the Ei(r'i) are real. If interaction is neglected, the trans- 
mission between the apertures is readily found 712 from the results of 
Hu" and Kay 5 : 



1 I f 1 f 1 
r - i / / Fadridn 

*-> \ Jo Jo 



(2) 



where 

Fu = -Si(ri)^2(r 2 )J"o(wir 2 )rir2 (3) 

and 

D = i- 2 T j l I Ei(n) I hidn £ \ E^) \ V 2 dr 2 l • (4) 

In these expressions the r, are normalized so that r, = rj/a,. The 
Fresnel number n = kaia 2 /d, with k the wave number, and J is the 
Bessel function of order zero. In the special case when the aperture 
separation is much greater than the aperture diameters, we see that 
n«l and that the aperture phases are uniform, i.e., &\{ri) = Ei{r % ). 
Substituting the small argument approximation for the Bessel function, 
x <£ 1, Jo(x) P3 1, eq. (2) then reduces to the familiar Friis trans- 
mission formula 16 

r-^f, »«i. (5) 

The effective aperture areas A\ are defined by 

f l &\{r)rdr * 
Jo 



A\ = 2ira 



/ $ 
Jo 



(r) 1 2 rdr 



i = 1, 2, (6) 



e.g., in the case of uniform illumination, 8\{r) = 1 and A\ = 7ra 2 <. 
The far-field transmission, T in (5), is also expressible in terms of the 
gains (G) of the apertures Ai and A 2 : 

T = G > G >(ird) 2 > n<<1 ' < 7 > 
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where 



0*-^r. i = 1 > 2 - W 



Returning now to the discussion of phase errors, suppose that the 
phases in the apertures depart from the ideal (confocal) distributions 
by amounts tfn(ri) in A i and <£ 2 (r 2 ) in A 2 . The transmission T n between 
the apertures in the presence of these phase errors is, from (2), 



1 I f l f 1 
Tu = j;\ / Fu exp lj(<f>i + 2 )]dridr 5 

Lf \ Jo Jo 



(9) 



where the &(r<) are abbreviated to <f>i. T u is expressible as 

Tu = To- A!T 12 , (10) 

where T is the transmission in the absence of phase errors and ATi 2 
is the loss resulting from the phase errors. Let 

Ti = To- AT 7 ,-, i =- 1, 2 (11) 

be the transmission between the apertures when there is a phase error 
on only the aperture At. From Appendix A we then find 

AT 12 ^ ATi + AT 2 + R, (12) 

where 

A7\ = I | J 1 J 1 Fudndn C f Fntfdndri 

- { f o f o F^idndrtY] , i=l, 2, (13) 

R = j: \ / Fi 2 dridr 2 / / Fufafadridr-i 

- f f Fi2<t>idridr 2 J J Fi 2 2 dridr 2 • (14) 

Equation (12) states that the total loss incurred by (small) phase 
errors <£i(ri) and 2 (r 2 ) on confocal apertures A i, As is approximately 
equal to the sum of the losses associated with each aperture when the 
other is free of phase errors together with the term R. In the next 
section we evaluate the expressions for A7\ and R when the aperture 
amplitude distributions are gaussian. 

III. GAUSSIAN APERTURES 

In the case of gaussian amplitude distributions, E\(?\) = exp (— a.Tj). 
To simplify the analysis we assume the a, to be sufficiently large that 
the upper limits in the integrals may be extended to infinity. The 
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transmission T in the absence of phase errors may now be obtained 
from (2), (3), and (4) by noting 16 

f°° 1 / n 2 r\ \ 

/ exp (— a 2 r 2 )J o(nrir 2 )r 2 dr 2 = ^- exp f — -j— J (15) 

to give 

T ° = 7JT7 v! » ex P (-«*) « x » * = X ' 2 - ( 16 ) 

(rr -f- 4ai«2) 

Apart from differences in notation, this expression is identical to the 
corresponding result obtained by Kogelnik 17 for the coupling of 
(fundamental) gaussian modes. We find 

T = 1 when n = 2^a^%, (17) 

i.e., within the approximation exp (—on) <K 1, there are optimum 
amplitude distributions that will ensure complete power transfer 
between the apertures for a given n. A detailed analysis, 9 or numerical 
integration, indicates this to be a satisfactory approximation when 
a» ^ 2.3, i = 1, 2. For example, when ai = a 2 = 2.36, the exact 9 
results for identical apertures are T = 0.9931, n = 5.00, and the 
approximate results are T = 1.00, n = 4.72. From (6), the effective 
aperture area, A% of a gaussian aperture is 

A? = — f , z=l,2. (18) 

Oti 

As expected from physical considerations A\ decreases as a, increases, 
i.e., as the aperture field becomes more concentrated about the aperture 
center. 

In the case of transmission with circularly symmetric, periodic phase 
errors, we take 

0<(r<) = /3,cos (jiTi), i = 1, 2, 
where (19) 

2irdi 



fii = k8 it ji = 



U 



Pi is the peak value of the error in radians (with 5, the peak profile 
deviation) and U is the period of the error. For errors of period much 
greater than the circumference of the apertures, 7, « 1 and then, in 
(9), 0»(r»-) £3/3,-, i = 1, 2. It follows, therefore, that T i2 = T , i.e., 
to this order of approximation, small, slowly varying, circularly 
symmetric phase errors do not affect transmission between the aper- 
tures. In the general case of small errors, the transmission loss AT\ 2 
for gaussian apertures is found from (12) with (13) and (14) by 
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substituting for the -E\(r,) and 0;(?',)- From Appendix B, we have 



^ = ^;[2D( 7 ;/2) - D( 7 D - 7^(7^2)], * = 1, 2, (20) 
i o 



where 



'* = * Vt 



Act, 



(21) 



'n 2 + 4aia2 
and 

S)(x) = exp (-r 2 ) Texp (r 2 )dr (22) 

is the (tabulated) Dawson integral. 18 The index y = {?} when i = {I}. 
The term R in (14) may be evaluated approximately in two cases 
of practical interest. In the first of these, the apertures are sufficiently 
far apart that n « 1 so that J^{nr\r%) ^ 1 in (3). Fu is then separable 
in functions of r\ and r 2 and hence, from (14), R = 0. For this case, 

AT 12 « AT 7 ! + AT 2) (23) 

i.e., the total loss is approximately the sum of the losses associated 
with each aperture when the other aperture is free of phase errors. 
The total loss is given by (20) and (23) in which y\ simplifies to 

tI«-&, n«l. (24) 

It is noted from (7) and (11) that 

^-75-', n«l, i=l,2, (25) 

where G\ = Gi — AG< is the gain of aperture A, with the phase error 
<f>i. Hence, (20) with (24) gives the fractional reduction in gain of the 
aperture A { resulting from a (small) periodic phase error. 

The second case of practical interest arises when the amplitude 
distributions on the apertures are optimized in accordance with (17) 
such that, in the absence of phase errors, the transmission is unity. 
From Appendix C, th e term R in (12) is negligible in this case provided 
7i, 72 2> n = 2Vaia2, i.e., the periods (Z,) of the phase errors satisfy 

h«- and h«— • (26) 

The transmission loss, ATi 2 , between the apertures is then the sum of 
the losses associated with each aperture as given by (23). This result 
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implies that the transmission through a sequence of confocal lenses, 
each with small phase errors of period satisfying (26), may be obtained 
by calculating the transmissions associated with each lens in the absence 
of phase errors on the other lenses. Furthermore, (23) indicates that 
it is not possible to compensate for such phase errors on one lens by 
introducing phase variations on an adjacent lens. When (26) is satisfied, 
the Dawson integrals in (20) may be replaced by the first terms of the 
asymptotic expansion (44) to give 

A7\ « W% i = 1, 2. (27) 

It is of interest to note the physical significance of the condition 
(26) for the validity of the approximate forms (23) and (27). As 
expected from the theory of diffraction gratings, a circularly symmetric, 
periodic phase perturbation on a circular aperture generates 19 two 
additional side lobes in the aperture radiation pattern. If the period of 
the phase perturbation is I, then the two side lobes are symmetrically 
located about the main beam at an angle 6 = sin -1 (X/l). Consider 
now the two apertures of Fig. 1 with phase errors of period h in 
Ai and h in A 2 . If the apertures are sufficiently far apart the side 
lobes, due to the phase error h in A\, will not intercept A 2 provided 
sin -1 (\/h) » tan -1 (a 2 /d), i.e., for small angles, h « Ad/a 2 . Similarly, 
the main beam of A\ will not couple energy to the side lobes of A 2 
provided Z 2 « \d/a\. The condition (26) implies, therefore, that energy 
is coupled from Ai to A 2 via the main beams alone. 

Figure 2 shows the transmission, as a function of y = 71 = 7 2 , 
between two identical apertures as obtained by numerical integration 
of (9) with (19), and as obtained from the approximate result (23) 
with (27). The upper curve in the figure applies for a = 4, n = 8, 
= 0.36 and the lower curve for a = 2.36, n = 5, = 0.18. In the 
absence of phase errors T = 1 for these (optimum) distributions. The 
dashed lines correspond to the. approximation (23) with (27). As 
anticipated earlier, the transmission is essentially unaffected by phase 
errors of large period, e.g., when 7 < n/2, i.e., I > 2Xd/a. The approxi- 
mate form (23) with (27) is seen to be within about 1 percent of the 
exact result when 7 > 2n, i.e., I < \d/2a. As an illustrative example, 
consider a beam waveguide system of the type described by Arnaud 
and Ruscio 2 with X = 3 X 10~ 3 m, d tt 80m, a ^ 0.5m. The parameters 
of this system correspond to those of the lower curve in Fig. 2. Substi- 
tution shows that small, circularly symmetric phase errors of period 
/ > 2a on the lenses will cause negligible transmission loss, and that 
the approximation (23) with (27) is applicable for phase errors of 
period I < a/2. 
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Fig. 2 — Dependence on y of phase error loss. 

IV. COMPARISON WITH PREVIOUS RESULTS 

To conclude this discussion on the effects of phase errors, we briefly 
compare the preceding results with the work of others. An expression 
for the gain ((?') of an aperture with small, periodic phase errors was 
given in (25). Consider the special case in which the period of the phase 
error is much less than the dimensions of the effective aperture, i.e., 
I « '4A*. From (18) and (19), this implies y » Va so that the Dawson 
integrals in (20) with (24) may be replaced by the large argument form 
(44) to give, with (25), 

% « 1 - tf», (28) 

where G is the gain in the absence of phase errors. Since this result 
depends only upon the magnitude /S of the phase error, it is anticipated 
that it may apply to random phase errors with correlation lengths 
that are small compared with the dimensions of the effective aperture. 
Ruze 11 has examined the reduction in aperture gain caused by such 
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random errors and we compare (28) with his results. In the particular 
case of a sinusoidal surface error of rms value e on a parabolic reflector 
antenna, we have = k8 = 2vlfce. From (28), the gain with this small 
phase error is 

f—[-( 4 f)l- 

This expression, derived here for a sinusoidal phase error, is identical 
to that obtained by Ruze in the case of a random error. As noted in 
Section I, Yoneyama and Nishida 13 have examined the effect of random 
phase errors on lenses in a two-dimensional, confocal, beam waveguide 
system. Their approach is based on the concept of a statistical beam 
mode and this leads to a description of the field distribution, and 
transmission loss, in terms of an integral equation. A computer was 
used to solve the integral equation by numerical iteration from the 
solution in the absence of phase errors. It is interesting to find that 
the conclusions of their study, of a two-dimensional system with 
random errors, are similar to those obtained here for transmission 
between circular apertures with periodic phase errors. In particular, 
it was found that the transmission was not appreciably affected by 
phase errors with large correlation lengths and that the loss for a 
given error tended to a constant value for increasingly small corre- 
lation lengths. 

V. CONCLUSIONS 

We have examined the effect of small, periodic, radial phase errors 
upon transmission between two coaxial, circularly symmetric apertures 
with confocal phase distributions. Two cases of practical interest have 
been considered when the amplitude distributions on the apertures 
are gaussian. In the first of these the apertures are widely separated 
with phase errors of arbitrary period. The total loss is then the sum of 
the losses associated with each aperture and is given in terms of 
tabulated Dawson integrals. This result reduces to a known form 
when the periods of the phase errors are sufficiently small. The second 
case of interest applies to transmission through a beam waveguide 
system with imperfect lenses. When the periods (U) of the phase errors 
on the apertures satisfy U < a,/2, (i = 1, 2), where a, is the aperture 
radius, the total loss resulting from phase errors is approximately 
\(fi\ + /3i), where /3i, /3 2 are the peak phase errors in radians on the 
two apertures. A comparison, based on numerical integration, shows 
this to be within about 1 percent of the exact result in a typical case. 
Phase errors with periods /, > 2a, (i = 1, 2) have comparatively little 
effect upon transmission. 

SMALL-PHASE-ERROR EFFECT 791 



VI. ACKNOWLEDGMENTS 

The author thanks S. J. Buchsbaum for suggesting this problem. 
Discussions with J. A. Arnaud and D. Gloge are greatly appreciated. 
The numerical integrations were kindly computed by D. Vitello. 

APPENDIX A 
Derivation of (12) 

For small phase errors, (0i + 02) <K 1 and the exponential in (9) 
may be expanded to second order. Recalling that the Ei are real, i.e., 
Fu is real, we then find 

t u « j) [ I SI f Fl2 ^ " * (01 + <t,i)2 ^ dTldr2 \ 

+ I £ f X F u {4>i + 0a)<Mr 9 rl ■ (30) 
Expanding the first bracket and noting that 

h\f f Fn< " <l>1 + <t>2 ^ dridr2 1 - ( * 1 + +*)$»» T °> < 31 ) 

where (</>i + 02)max is the maximum value of (0i + #2), we have, to 

second order in 0i, 02, 

T 12 = T - ATn, (32) 

where 

ATi2 » 5 [ f f Pndndn C C F u (4>i + to)Hr x dn 

- j £ £ Fn(4n + H)dndrX j • (33) 
Expanding the brackets gives (12) with (13) and (14). 

APPENDIX B 
Derivation of (20) 

Substituting £<(tv) = exp (— a.-r?), 0,(r.) = frcos (yin) (»' = 1, 2) 
into (13) with (3) and (4) gives, with (15), 

^ = ariCIiM - 2„/f( 7l )] i - 1, 2, (34) 

where 

h(yd = / exp ( — rjr 2 ) cos 2 (yrfrdr, (35) 

./o 

JaM = / exp ( — ?jr 2 ) cos (jir)rdr (36) 

Jo 
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and 

v = £ (n 2 + WO, (37) 

withy = {?} when i = \\}. Expanding cos 2 (7,7-) and integrating: 

hM = 5 + |/ 2 (27.). (38) 

Integrating by parts, 

h(y<) - ^ Fl - 7. | o " exp (-r,r 2 ) sin (7<r)drl , (39) 

which is expressible 18 in terms of the (tabulated) Dawson integral, i.e., 

'■w-il 1 -$*(&)]■ <40) 

where 

»(x) = exp (-a; 2 ) f" exp (T 2 )dr. (41) 

From (34), (38), and (40), we then obtain (20) and from (19) and (37) 
we obtain (21). 

APPENDIX C 

Approximate Evaluation of R In (14) 

We derive an approximate expression for R when 71, 72 ^> n. It is 
assumed that the amplitude distrib ution s on the apertures are opti- 
mized such that T = 1 with n = 2Vaia 2 , where oi, 0:2 ^ 2.3. Consider 
the integrals in (14) : Extending the integration limits to infinity and 
substituting Ei(n) = exp (— a,r, 2 ), * = 1, 2 gives, with (3), (4), 
and (15), 

/.' /„' *»*•*' - h • D = h- < 42 > 

Similarly, substituting </>, = /3, cos (7,7-,) and using (15) and (40), 

/;/;^ lWr ,-A[ 1 _i^ D (^)]. m 

Since 71 )$> n, the Dawson integral may be replaced by the first two 
terms of the asymptotic expansion, 20 



to give 



f f Futpidridr* « - jr^-g- (45) 

'0 7o ^0:271 
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Further, 

/ / Fnfafadridrt = /3i/3 2 / exp (-airf)ri cos (yiri)3(ri)dr h (46) 
Jo Jo Jo 

where 

0(n) — / exp (— a2HD«/o(nrir 2 )r2 cos (y^dn. (47) 

Substituting the integral representation of the Bessel function 

1 f 
J (x) = - cos (a; sin d)dd, (48) 

IT Jo 

interchanging orders of integration and expanding the cosine product, 

tin) = ^ f* [/(*) + K-e)2de, (49) 

where 

1(6) = / exp (— a 2 r|)r 2 cos [(72 + nr\ sin 0)r 2 ]dr 2 . (50) 

From (40), 

*(») = JT- F 1 - "7= (Y8 + «ri sin 0) 
■6«2 L Va 2 



° 2 

—7= (72 + fifi sin 6)\\> (51) 



Since 72 » n, both Dawson integrals in (49) may be replaced by the 
large argument form (44) to give 

din) tt 1(d) tt-y* 2 . (52) 

Evaluating (46) by (40) and using (44), 

f [ Fufh^ridrtPsi^ri- (53) 

Jo Jo 7i72 



Substituting for the integrals in (14) and reducing (20) then gives 
R _ ftp, 8n 2 



(54) 



A7\ + AT 2 « + «7!7i 

But |/3ij8,/(/8? + /3S)| ^1, i.e., 

I / On \2 

(55) 



fl l^/2n_V. 



ATi + AT 2 | - V7172 
Since 7i7 2 2s> 2n, (55) is much less than unity and so, from (12), 

AT 12 W AT 7 ! + A7V (56) 
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